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Abstract. We examine the possibility of dealing with gravitational singularities 
on a quantum level through the use of coherent state or wavelet quantization in- 
stead of canonical quantization. We consider the Robertson- Walker metric coupled 
to a perfect fluid. It is the simplest model of a gravitational collapse and the results 
obtained here may serve as a useful starting point for more complex investigations 
in future. The coherent states, or equivalently the continuous wavelet basis ele- 
ments, are built from the unitary irreducible representation of the affine group of 
the real line with positive dilation. A natural regularization of the singularity is 
an interesting issue of our approach. 
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1. Introduction 



The purpose of this article is to examine the implementing of the coherent state 
quantization in the study of the gravitational singularities. We consider the Robertson- 
Walker metric coupled to a perfect fluid. It is the simplest model of a gravitational 
collapse and the results obtained here may be a useful starting point for more com- 
plex investigations in future. 

Canonical quantization of a Friedman-Lemaitre universe with an eye towards the 
fate of gravitational singularity has been studied extensively. Early treatments in- 
clude Blyth and Isham [T], featuring a discussion on the ambiguous meaning of 
singularity resolution, and the work [2] by Lapchinskii and Rubakov, who obtained 
a quantum non-singular perfect fluid-filled universe. However, the interpretation 
of these and other, more recent, results is not obvious for at least three reasons. 
Firstly, in [2] it was shown that the classical and singular evolution can be replaced 
by a unitary and thus non-singular one, provided one fixes an appropriate boundary 
condition to ensure self-adjointness of the Hamiltonian. Unfortunately, there are 
infinitely many 'equally good' choices for the boundary condition. The choice has to 
be made without a clear justification. Secondly, at the fundamental level, the pre- 
vailing attitude towards quantum gravity is that one needs some additional input, 
let it be loops, strings or triangulations. The inclusion of these basic premises should 
give rise to novel effects in the context of mini-superspaces. Finally, there is the 
so-called 'problem of time': in canonical quantum gravity one is forced to describe 
the evolution of the space with respect to a chosen degree of freedom and it may be 
shown that the resultant quantum theory significantly depends on this choice [3] . In 
this light, the expectation that Weyl quantization of mini-superspaces should lead 
to some important clues about the singularity may turn naive. 

Recently, it has been argued that quantum cosmology to become an empirical 
science can only be developed as an effective theory [3]. There are simply too many 
unknowns, including unresolved conceptual issues and technical complexities, to have 
a hope that a fundamental and rigorous approach is feasible. Bojowald suggests that 
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a good cosmological theory should therefore be flexible enough to parameterize our 
ignorance. Canonical (i.e. Weyl) quantization is in a sense rigid and will not provide 
us with parametrizable physical models. For a clear and comprehensive review of 
various quantization methods (e.g., canonical, geometric, deformation, ...) we refer 
to Ali and Englis in [6J. In what follows, we propose to relax the usual, canonical, 
quantization prescription by implementing coherent state quantization in the study 
of the cosmological singularity. 

The coherent state quantization was demonstrated to be a valid alternative to 
canonical quantization in dealing with various simple systems (see (3, |8] and refer- 
ences therein). It is a flexible method, because it allows for a reasonable amount of 
freedom due to the free choice of coherent states (or 'wavelet basis') which determine 
the quantum realm of a model. In CS quantization ([§], see also PJ for deep proba- 
bilistic aspects of the procedure and for a various list of elementary examples), one 
does not require that the Poisson bracket of basic variables is strictly mapped into 
the corresponding Lie algebra of the operators. Besides the standard linearity, iden- 
tity correspondence 1 i— > I and the self-adjointness of quantum observables, the only 
minimal conditions are: (i) the quantization must agree with available measurements 
and (ii) it should consistently admit the classical limit. In this way one allows for 
many more possibilities. This seems to be a desired quality while quantizing gravita- 
tional fields. Surprisingly, relaxing some of the usual quantization constraints leads 
to the occurrence of a quantum repulsive potential, which regularizes the singularity 
and leads to a unique unitary evolution across the big bounce. This is the central re- 
sult of this work. Apart from that, a quantum 'true' Hamiltonian, which is obtained 
by the intuitive construction, is naturally associated with a semiclassical counterpart 
through the so-called 'lower symbol'. Last but not least, CS quantization map has 
the advantage of being well-suited to deal with non-standard (e.g. non-polynomial) 
Hamiltonians occurring in gravitational systems, which may be useful in future in- 
vestigations. 

This paper is organized as follows. In Section 2 we recall briefly the equations of 
the FLRW model, we define the Hamiltonian and the relevant physical quantities. 
In Section 3 we present the details of our quantization scheme based on an Unitary 
Irreducible Representation (UIR) of the affine group. At this stage some important 
parameters of the representation are not specified. The expressions of quantized 
observables are given in Section 4. In Section 5 the parameters of the representation 
are specified, we obtain analytical expressions, and we perform numerical simulations. 
We analyze the physical aspects of our results in Section 6 and we conclude in Section 
7. The appendix is a short review of the affine group and its representation (s). 
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2. The FLRW model 

2.1. Hamiltonian constraint. The space + time split in the Hilbert-Einstein ac- 
tion leads to the Hamiltonian density as a sum of first-class constraints [TU] : 

(1) H = NC° + NiC 1 , 
where 

(2) C° = -q x ' 2 ( 3 R + q- ll ] 



Generalized positions q^ correspond to the three-metric on spatial sections of space- 
time, 3 R is the three-metric Ricci scalar, and p l i are generalized momenta, associated 
with extrinsic curvature of the three-hypersurfaces. The space {(qij,P^)} is equipped 
with the Poisson structure 

(3) {q ij (x),p kl (x')} = 8 { *5 j \8(x-x'). 

Scalar N and covector iVj encode the remaining components of space-time metric, 
which here play a role of Lagrange multipliers. In the above we have put 

4 k '.— — =— « 1.24 x l(T 34 skg _1 = 1 . 

c 6 

The units will be restored if necessary, within the quantum framework when we 
substitute H y 1/kH while inserting the Planck area-like dp = kH ~ 8.2 x 10 _68 m 2 
instead of h. This a p yields the natural length standard v/ap. 
In the FLRW universe the line element is given by: 

(5) ds 2 = -N{t) 2 dt 2 + a(t)%wV , 

where 1-forms id 1 are invariant with respect to the homogeneity group of spatial 
section. The specific form of u l depends on the spatial curvature (see e.g. [H]). For 
metric tensor's components in ^ we introduce: 



/ adu _ -12 



(6) a := — -y, p a := — -y- / a—du, {a,p a } = 1, a > 



(Jdco)i- (Jdu) 



N 



where du = u 1 A u 2 A u 3 . The integration is performed over the whole universe if 
it is compact and over any finite patch otherwise. Thanks to the definition of basic 
variables ^ we get rid of the Dirac delta featuring in (|3j), identify the three non-zero 
metric components, give a the interpretation of a distance and make it convenient 
to express (the integral form of) the Hamiltonian in terms of them. 

The vacuum formulation ([lj [2J was extended to include perfect fluids by Schutz 
[12] . He used the potential- velocity formulation: u v = j^(4>,u + ot^ tV + 9s }U ), where 
u v are the components of the proper velocity of the fluid's element, (0, a, (3, 9, s) are 
independent potentials and h is the specific enthalpy, which is not an independent 
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variable due to u u u v = — 1. The specification of the indicated framework to the 
FLRW models filled with barotropic fluid subject to the equation of state p = wp 
(w = const) is straightforward and may be found in [2J. In this case the Hamiltonian 
constraint is 



(7) H — N (-^ - 6ka + « . 
w V 24a a 3w J 

where the dimensionless k = (J dco) 2 ^ 3 k and k = 0, — 1 or 1 depending on whether 
the universe is flat, open or closecQ Parameters (T,px) G 1R x M.+ are a canonical pair 
associated with the fluid. More specifically, the meaning of pt is the fluid's energy 
times the volume in which it is contained to the power of w and T is an auxiliary 
function of the physical dimension of length to the power of 1 — 3w. The dimensions 
will explicitly agree in (J7) once we restore - in front of the gravitational part of the 
Hamiltonian constraint, n order to bring (ffl) to more convenient form we define new 
basic variables for geometric observables: 

(8) (q, P ) := (a^'\ -^^fi^-w) , {q,p} = 1, q > . 

Note that the physical dimensions of q and p: 

(9) [g] = L 3 ( 1 -^ 2 , [p] = L (3u,+1 )/ 2 , 

so that [qp] = L 2 . After specifying the lapse iV = g 2u, /( 1 ~ u ') ) we eventually obtain: 
(10) 



H = (-a(w)p 2 - 6kq^ w) + Pt ) w , p:(w) 
where [H] = L^ 3w+1 \ 



2(3w + l) . , 3(1 -w) 
a[w) " 



3(1 - w) y ' 32 



2.2. Reduced phase space. The Hamiltonian formulation of general relativistic 
models introduces first-class constraints, given in (J2j), which reflect the coordinate 
freedom in the Einstein theory. Here, the physical symmetry of the FLRW models 
enables one to make use of the preferred foliation of space-time and reduce the 



formulation considerably to a single constraint (10). To implement quantization one 
may employ either the Dirac approach ('first quantize, then solve constraints') or 
the reduced phase space approach ('first solve constraints, then quantize'). However, 
both the constrained and the reduced phase space will accommodate an incomplete 



1 We note that the value of J dco is not arbitrary for 'a' has the geometrical meaning of the 
radius of curvature. For example in the case of the 3-sphere (k=l) we have J du = 2ir 2 . 
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dynamical flow due to the singularity (in a suitable choice of time variable)^] In what 
follows we restrict ourselves to the reduced (unconstrained) phase space analysis, with 
the choice of the variable 'T' as a clock. 

The reduction of the model goes this way: we start with reducing the symplectic 
form O to a closed two- form Q^: 

(11) tt = dqAdp + dT A dp T ->■ n R = dqAdp + dT Ad (a(w)p 2 + 6kq» (w) ^ . 

where we solved the constraint for p?. The form Qp lives on the constraint surface 
and is not symplectic, because it is degenerate with the Hamiltonian vector field 
Vh = {•, H} being in its null direction. In order to get the physical Poisson bracket 
we put T = const and then invert Qr. Now our system, reduced and no longer 
constrained, is given by: 

(12) {q,p} = l, h T = a(w)p 2 + 6kq^ w) , q>0. 

Hence we model a singular universe as a particle moving on the half-line, where 
the end-point of the half-line signals singularity at which the classical dynamics 
terminates. From the explicit form of the Hamiltonian it is apparent that the clock T 
is slow-gauge^J Two cases are particularly interesting: (i) w — 1/3, which corresponds 
to the radiation as the content of universe; (ii) k — 0, which corresponds to flat FLRW 
universe, which apparently can be modeled as a freely moving particle. 

The compound geometric observables of physical interest include the volume V 
and the expansion rate 6, i.e., the trace of extrinsic curvature: 



1—w 



(13) V:= / a 6 dio = q^, 6:=- = -(1 - w)pq 

J 1 — w IN q 8 

As the singularity is approached V — > and 9 — » ±oo. 

3. Quantization of the half-plane 

Because we are going to use affine transformations of physical quantities, we should 
keep control of the physical dimensions. In view of this, we introduce the parame- 
ter u = ap 1 which has the physical dimension of q. The dimensionless scale- 
momentum half-plane is then defined as n + = {(q, p) | p e K , q > 0}, where 

(14) q:=^ p:=*±- 

a ap 



2 One may argue that even for non-singular models there may exist a choice of clock with an 
incomplete Hamiltonian flow. Here we exclude such clocks. For a discussion of the dependence of 
quantum theory on the choice of time (so-called 'multiple choice problem') we refer reader to [3]. 

3 A clock is slow-gauge if the singularity is reached within a finite time interval [5]. 
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Equipped with the multiplication 

(15) (q, p) • (q , Po) = (qqo, - + p) , q G K, p e R , 

it is viewed as the affine group Aff+(R) of the real line (see Appendix for more de- 
tails). This group possesses two non-equivalent unitary irreducible representations 
(UIR) U±, besides the trivial one. Both are square integrable and this property 
is fundamental for the continuous wavelet analysis [HI [151 HSJ E3 HI US!- Equiv- 
alent realizations of one of them, say U + = U, are carried on by Hilbert spaces 
H a = L 2 (R* + ,dx/x a+1 ) where a is an arbitrary real number. Nonetheless, detailed 
calculations prove that these multiple possibilities do not introduce interesting dif- 
ferences. Therefore we choose in the sequel the standard case a = — 1 and denote 
U = %-x = L 2 (lR;,dx). The UIR of Aff+(R), expressed in terms of the physical 
phase space variables, acts on % as 

(16) U(q,p)^(x) = ^ x ^(x/q). 

vq 

Given a normalized vector ipo e "H, a continuous family of unit vectors are defined 

as 

(17) \q,p) = U(q,p)\if> ), (x\q,p) = e ipx ^-^{x / (\) . 

vq 

where orthonormal basis \x) in a distributional sense obeys (x\y) = 5(x — y). The 
transported ipo is the (mother) wavelet or fiducial vector. For the sake of simplicity 
we assume in the sequel ipo(x) is real- valued. 
Let us define the constants c a (for real a) as 

(18) c a := [ ip (x n 



x 2+a ' 



the resolution of the unity is straightforward: 

f dqdp 1 

(19) / \q,p){q,p\ = l, 



provided that c_j < oo. Due to this crucial property, the family (17) is called 



continuous wavelet basis (in signal analysis) or coherent state family for the affine 



group within a more quantum oriented context. Thanks to (19), the CS quantization 
of classical functions f(q,p) can be implemented through 

(20) f*A,= I ^^f(q,p)\q,p)(q,p\. 
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Indeed, when properly defined, this map is linear, to / = 1 there corresponds the 
identity operator and to real semi-bounded / there corresponds self-adjoint operator, 
i.e. three basic requirements of any quantization procedure. 

4. Physical operators 
4.1. Functions of q. The quantization of coordinate functions reads as 

(21) q ^ A q = —Q , Q<f>(x) := ax<f>(x) , 

provided that c < oo. This operator is self-adjoint. 
More generally for e M 

(22) q P^ A c lz±Q^ 

C-l 

provided that cg_x < oo. 

2 

In particular the quantization of the volume V = gi-«> reads as 



(23) A v = ^^Qi-™ 



C l+w 2 

1 — w 

C-l 

4.2. Functions of p. The quantization of momentum reads as 

(24) p^A p = P, P(f>(x) := -i—4>'(x) 

a 

This operator is symmetric when defined on rapidly decreasing functions with 
support in R^, but one can show that there is no self-adjoint extension (the deficiency 
indices are (1,0) [20]). 

For the "/ree" Hamiltonian Hq = A, p i we have: 

K f°° vr\v 

(25) Ho = P 2 + a 2 P —, with K = K(ik):= {^{u)) 2 

Q 2 Jo c -i 

We notice that the quantization procedure always yields an additional term. This 
term depends only on the fiducial vector and its importance will be explained below. 



4.3. The quantized Hamiltonian. The classical Hamiltonian of (12) reads as 

(26) h T (q,p) = a(w)p 2 + 6kq^ w) , with jfe = (/ dw) 2/3 jfe, k = 0,±1 , 
therefore the quantized Hamiltonian H = Ah T reads as 

(27) H = a(w)P 2 + a 2 P a(w)^ + Gk^^Q^ . 

Q 2 c_i 

Let us notice that K > 0, therefore the singularity x = cannot be reached: this 
dressing of the classical singularity is an outcome of the followed CS quantization 
scheme (through K). Furthermore our procedure yields a renormalization of the 
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coupling constant of the potential. 

If we assume a closed universe with a radiation content then w = 1/3, k = +1 and 
k = (J dw) 2//3 , then variables g and p get both a length dimension (Eq.([9j)) and we 
obtain the special Hamiltonian H cr 

(28) H , = ^ + ^ +ei ^. 

This Hamiltonian is an ordinary differential Sturm-Liouville operator, singular at 
the end point x = 0. The functional properties depend on the value of K, as follows 
from the analysis of Gesztesy et al [21] (see also [2D])- In particular, K = 3/4 is 
the critical value, while one would naively expect K = 0, i.e., the infinite barrier, to 
play the role. Using the standard approach and terminology of [2D] , the potential 
term in H cr is in the limit point case at the end point x = 0iffT>3/4 and in the 
limit circle case atx = 0if0<-ft'<3/4. The potential is in the limit point case at 
infinity. It follows that H cr (defined on the domain of smooth compactly supported 
functions) is essentially self-adjoint in the former case. In the latter range of K, the 
deficiency indices of H cr are (2, 2) and therefore more self-adjoint extensions exist; 
see [21] for a detailed analysis. In this paper, we choose the fiducial vector ipo in 
such a way that K > 3/4 in order to comply with essential self-adjointness of the 
Hamiltonian. This ensures non-ambiguous time evolution and singularity resolution 
at the quantum level. 

4.4. The quantized expansion rate O. We have 

Q 3C 2_ 1 1 

(29) 6 = -(l- w)pq-T^ ^ Ae = Q = -^(1 - w)^^P^^. 

O O C_l Q2(1-to) Q2(1-«i) 

This operator is symmetric when defined on a suitable domain, but it is not self- 
adjoint (or does not possess a self-adjoint extension), from the reasoning previously 
developed for P. 

5. Explicit formulae in a particular case 

In the previous formula the wavelet ipo is a free parameter of our CS quantization. 
In order to obtain explicit expressions, we decide to fix ip as the following unit vector 
in^: 

(30) ipoH%) = 1 e ~^ x+ ~k), with v > and f > 0. 

y/2xK (v) 

We notice that i/jq (x) falls off with all its derivatives at the origin and at the infinity. 
To calculate the normalization constant involved in the previous definition of i/j^ 
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(and in order to obtain other useful integrals) we take benefit of the formula [22] 

(31) Va, b, c e C, n{b) > 0, sR(c) > 0, J x^e'^^dx = 2 I- K a {2\/bc) . 

Here, K a is a modified Bessel function |23j . We recall that its asymptotic behavior at 
large argument v is K a {v) ~ e~ v ' 'it / '(2is) , whereas at small v <C \J a + 1, K a {y) ~ 

(l/2)r(a)(2/z/) a for a > and i*T (f) ln(z//2) -7. 

The coefficients c Q defined in Eq.(18) reads as 

(32) c a = p^-K a+2 (u). 
The coefficient £ is fixed in the sequel as being 

(33) e - 

The interest of this choice will become apparent with the study of the operator A q . 
Also note that it implies that, in view of the study of the semi-classical regime, ap 
and v remain the only free parameters. 
The coefficient c_i involved into the resolution of unity is 

[ } 1 ~ K (u)K 2 (uY 

5.1. The physical operators. The quantization of the variable q reads as 

(35) A q =^Q = Q, 

due to of the choice of £. The quantization of the potential q@ leads to 

(36) A q , - c i Q - K ^ Q . 

This formula can be applied in particular to the classical volume V = q 2 /( 1 ~ w \ 
The quantization of the momentum P = A p is as in (24) while the "free" Hamiltonian 
H = A p 2 becomes 

(37) ir = p. + 4*M, AT M :=ifl + „*■<"> 



The Hamiltonian if of Eq.(27) reads as 
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The last quantized observable is the expansion rate 6 = Ag that reads as 



l+w 



e 



s 



l+w 



-p- 



l+w 



(39) _ 

K X (p)^ Q2(l-«) Q2(l-v>) 

Finally, a special interesting case of Hamiltonian H is obtained when we assume a 
closed universe with a radiation content (k = +1, w — 1/3 and k = (J dcu) 2//3 ). Then 
the Hamiltonian H cr of Eq. ( 38 ) reads as 



(40) H cr 
5.2. Analysis of H cr . 



24 



-P 2 



ap K(u) 

24 "g 2 " 



+ 6A; 



KM* 



Q 2 . 



5.2.1. Semi- classical point of view. If we perform a semi-classical analysis of the 



dynamics due to H cr obtained in Eq.(40), we see that the supplementary repulsive 



potential generated by the quantization leads to a displacement of the equilibrium 
point of the potential. While the harmonic potential alone possesses an equilibrium 
point located at the singularity q = 0, the harmonic potential corrected with the 
supplementary repulsive term exhibits a different equilibrium point which is located 
at 



(41) 



v) 



Qe 



K{v) 



We recover q e = when ap = 0, but for ap ^ we find that the smallest value 
lip/ 5.8 is obtained for v — > 0, while q e — > oo when v — > oo. Otherwise the 



q e = 

renormalized coupling constant of the potential in Eq.(40) converges to the classical 
expression when v — > oo. Therefore the free parameter v and the standard ap can 
be used to specify the renormalized coefficient of the potential and the position of 
the new equilibrium point. 

We observe that the semi-classical dynamical behavior of the system is an oscillation 



around this point q e and the "bare (or true) classical singularity q 
reached. 



V 



is never 



5.2.2. Eigenstates and evolution operator. With the quantities £, u and A defined as 
(42) i = 




kK x {v)Kz{v) 



/v,(>) " / ' " K 2 {v) 
the eigenvalues E n , n = 0, 1, . . . of H cr read as 



, and A = 6u . 



(43) 



En 



a P cu(2n + £ + 3/2). 
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The corresponding normalized eigenvectors ip n are 
(44) ij n (x) = N n x e+1 L e n +1 / 2 (2Xx 



2 )e- Xx 



The functions L^ 1 ^ 2 are the associated Laguerre polynomials [23] and the normal- 
ization factor N n is given by 

, 2(2A)^+ 3 / 2 n! xV2 
(45) N n -' 



T(n + £ + 3/2) 

Introducing a dimensionless evolution parameter r related to time t through some 
scaling for instance, the evolution operator U (r) = Q- lH crT/a F j g p er i oc ij c with period 
ir/u. It is given in terms of matrix elements (x\U(r)\y) by: 



(46) (x\U(r)\y) = 2A , ^^ e ^W)^) j I 



sm{uJT)\ \\sm{uT) 

where J v is a Bessel function. This expression is derived from series involving 
Laguerre polynomials [22J. The r.h.s. of Eq. (46) is singular for sin(wr) = 0, 
(x\U(r)\y) being in that case the distribution 5{x — y), corresponding to U(r) — 1. 

5.3. Lower symbols. The expectation values (q,p\A\q,p) of quantum operators al- 
low to map the quantum world into the classical one and they are called "lower sym- 
bols" [21] or "covariant symbols" [25J. Therefore we can examine the semi-classical 
map f(q,p) i-> f(q,p) = (q,p\Af\q,p) that exhibits the corrections (regularizations) 
induced by our quantization procedure. The map />->•/ reads as the integral trans- 
form 

(47) f{q,p)= [ ^L\( q ,p\ q >,p>)\*f( q >,p>). 

It is the average value of the function f(q,p) with respect to the probability dis- 
tribution (q',p') i — y 2napc i KQ'PW iP')\ 2 - Viewed as a kernel, the latter is expected 
to play a regularizing role under the form of a generalized convolution. Due to the 
resolution of the unity (19), the scalar product (q,p\q',p') is a reproducing kernel for 
a Hilbert subspace of L 2 ljl + , dqdp). It is given by 



/ km / i / n ■ r - / q + l' L , ^iqq'ip 1 -p)\ K^u) 
(48) (q,p\q ,p) = -^—^Kq v ^^\ 1 + <■( , a » Wlth £ " 



K Q {y) V SvWV apu^q + q'))' s K 2 (u) 
For the powers of q we have 

(49) ( q ,p\A„\ q ,p) = 
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in particular 

(50) (q,p\A q \q,p) = (q,p\Q\q,p) = ttj^q- 

We recover (q,p\A q \q,p) ~ q if v — \ oo. Otherwise we obtain for P = A p 

(51) (q,p\A p \q,p) = {q,p\P\q,p} = P- 
The expectation value (q,p\P 2 \q,p) reads as 

K x {v) 2 K{v) 1 



(52) (q,p\P 2 \q,p) =P 2 + a 2 P 



K Q (v)K 2 (v) q 2 



We notice the supplementary term in 1/q . Using the previous equations (36), (49) 
and (52) we obtain 

K x {v) 2 K{v) 1 



(53) (q,p\A p2 \q,p) = p + 2a P 



K (p)K 2 (p) q 2 



We notice that the 1/q 2 coefficient in (q,p\A p 2\q,p) and in (q,p\P 2 \q,p) differ from 
a factor 1/2. 

The coupling constant of the supplementary term in q~ 2 disappears if ap — > or 
v — 7- 0, it becomes infinite if v — > oo (for a fixed value of ap). Finally we obtain for 
the classical Hamiltonian of Eq.(12) 

(54) h T {q,p) = a(w)p + 2a P a(w) ? + 6fc K^K^u) 9 " 

We recover hx — > hp (and — > q 13 , p 2 — > p 2 ) if we assume first ap — > and 
then z/ — ?• oo (the limits are not commuting). If we want to keep two independent 
parameters and to have independent limits, we need to "renormalize" ap, including 
in its definition some of the i^-Bessel functions involved in the coupling coefficient 
of q~ 2 as for example 



a P K Q (v)K 2 (v) 

(55) ap 



2 V K x {v) 2 K{v) 



where ap would be the new renormalized Planck area. Another solution (to avoid 
the problem of non commuting limits) is to assume that v is in fact a function of ap. 
The expression viap) must be well chosen: we need to impose both a 2 P v[ap) — > 
and v{ap) — > oo when ap — > (a simple solution is v{ap) oc 1/ap). 



14 



H BERGERON, A DAPOR, J-P GAZEAU AND P MALKIEWICZ 



5.4. Time evolution in phase space. For any normalized state G "H, the resolu- 
tion of unity allows us to get its phase space representation $(g, p) = (q, p\<p) / yJ2-na P c_i 
and the resulting probability distribution on the phase space n + : 



(56) 



11+ 3 (q,p) h-> 



27rapc_i 



\(<i,p\0)\ 2 = p<p{q,p) 





IS 1-5 2.0 2.5 



Figure 1. Phase space distributions (56): on the left for the eigen- 
state ipo °f H cr , on the right for the coherent state \qo,po) with q Q = 2, 
Po = 0. The parameters are fixed to the values v = 80, ap = 1 and 
k — 1. The thick curve (figure on the right) is the expected phase tra- 



jectory, deduced from the semi-classical hamiltonian in Eq. (58). The 
ranges of variables q and p are respectively [0.2, 2.8] and [—35, +35]. 
Increasing values of the function are encoded by the colors from blue 
to red. 

In Figure [T] two phase-space distributions p<p(q,p) are shown: for the ground state 
ipo of H cr and for the coherent state \qo,Po) with (qo,po) = (2,0). In both cases the 
quantization parameters are chosen as v = 80, ap = 1 and k — 1. 

Let us now examine the time behavior r \-t p</,( r )(9,p) for a state <p(r) evolving 
under the action of the Hamiltonian H cr of Eq. ( 40 ) : 



(57) 



(r)) = e- lHcrT/ap \<p) 

oo 



n=0 



where \ ipn) is eigenstate of H cr with eigenvalue E n given in Eqs. (43), (44). 
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With |0) = | qo,Po) as an initial state and using Eq.(54), we have for a given v 
(58) (qo,Po,T\H cr \q ,p ,T) 



1 r? I n * MM 
24 Po + ap 12iro(^ 2 (^) % 2 



I , aZ K 2 {u)K 3 {u) 2 



Since for large values of v the lower symbols of Q = A q and P = A p correspond to 
their classical original functions q and p, one can expect that the time average of the 
probability law P\ qQ , P0 , T )(p, q) corresponds to some fuzzy extension in phase space of 
the classical trajectory corresponding to the time-independent Hamiltonian in the 



r.h.s. of Eq.(58). Similar dynamical issues of CS quantization are encountered in the 



case of Poschl- Teller potentials [7] . 




Figure 2. Time average phase space distribution of Eq.(59) for 
|9o;Po)( r ) evolving following the Hamiltonian H cr . The values of pa- 
rameters and the thick curve are those of Figjl] The ranges in q and 
p are respectively [0.2, 2.8] and [—35, +35]. The domain < q < 0.2 
is not represented because of numerical instabilities. Increasing values 
of the function are encoded by the colors from blue to red. 

This key result is illustrated in Figure [2] where we have represented the time average 
distribution p defined as 
(59) 

1 f T 1 °° 

p(q,p)=lim — I p\ Q0 ,p ,r)(q,p)dT=^——^2\(q,p\^ n )\ 2 \(^ n \q ,p )\ 2 , 



2napC-\ 



n=0 



for the same values of the parameters as in Figure [TJ The time average distribution 
p allows us to compare the quantum behavior with the classical trajectory, but the 



expression (59) hides the details of the wave-packet dynamics, i.e. the bouncing of 



the wave-packet during its periodic motion. The figure represents this behavior. 



16 



H BERGERON, A DAPOR, J-P GAZEAU AND P MALKIEWICZ 





Figure 3. Phase space distributions P\q , Po , T )(p, ( l) at different times 
equally spaced (from top left to bottom right). The values of parame- 
ters and the thick curve are those of FigjTJ The ranges in q and p are 
respectively [0.2, 2.8] and [—35, +35]. Increasing values of the function 
are encoded by the colors from blue to red. 



6. Discussion and Conclusion 

By the use of coherent state quantization we have derived a quantum model of the 
Friedman-Lemaitre universe, in which the classical big bang singularity is replaced 
with a quantum big bounce resulting in a smooth and complete evolution. At the 
most general level, the removal of singularity is due to a unitary evolution: a gravita- 
tional collapse, represented by the Dirac delta peaked at q = can never be reached 
as this state does not belong to the Hilbert space. The novelty introduced by CS 
quantization is that the singularity resolution is accompanied by the occurrence of a 
repulsive potential. The potential's role is twofold: (i) on the quantum level it may 
lead to a non-ambiguous unitary evolution across the bounce, (ii) on the semiclassical 
level it provides a mechanism for the universe to stop contracting, bounce and start 
expanding. The pre-factor in the potential term (which will be discussed below) is 
the desired free parameter that provides a way to match our theory to observational 
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constraints. In particular it may be used to set the energy scale at which quantum 
effects come into play. Let us emphasize that this inverse quadratic potential arises 
for all possible choices of the fiducial vector ip in our quantization scheme. 

6.1. Finite universe. In the case for which the universe is non-compact the Hamil- 
tonian formulation is derived through the restriction to a finite patch of space. The 
inspection of the semiclassical Hamiltonian and the definition of the basic variables 

and (|8) and the curvature constant k shows the following: both the kinetic and 
curvature terms of the Hamiltonian depend on the size of the patch like (J du) w+1 
while the potential term behaves like (J do;)™ -1 . In other words, the classical dy- 
namics is invariant with respect to the choice of the patch, whereas the quantum 
dynamics is not. The repulsive potential breaks the invariance and so the physical 
content of the theory depends on the patch. Thus, one has to exclude non-compact 
universe from quantum modeling. We note that it does not imply the curvature of 
the universe: for k = we may consider e.g. torus topology, for k = 1 the universe 
is necessarily compact, and for k = — 1 there are infinitely many possible compact 
spatial sections [26J. 

6.2. Planck era and the quantum phase. From the lower symbol of Hamilton- 



ian (54) we deduce that the quantum effects become important once the repulsive 
potential gets comparable with the kinetic part, that is p 2 ~ ((3 2 (v)/q 2 ), where 

/3» = 2a 2 P K 1 (u) 2 K(u)/(K (u)K 2 (u)) . 

In other words, the region of phase space in which classical dynamics cannot be 
trusted is given by the inequality (with curvature neglected): 

\pq\ < , 

which by making use of the definitions collected in the introductory section one can 
translate into the geometric observables: 

(60) \ve\ < . 

Thus the quantum effects do not depend neither on the size of the universe nor on 
its expansion rate alone but rather on the specific combination of them both. It is a 
common belief that the quantum dynamics begins roughly when the energy density 
of matter hits the Planck scale and the universe is said to enter the Planck era. 
However, our result does not confirm this assumption. First we note that by virtue 
of the Friedman equation the energy density of matter content is proportional to the 
expansion rate squared if the curvature term is neglected [^J However, inequality (60) 



4 This is a reasonable assumption, because as the universe approaches the singularity the intrinsic 
curvature term becomes dominated by the energy density of fluid. 
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says that even for the Planck scale value of expansion rate the universe may still 
be classical provided it is large enough. On the other hand, a low energy density 
universe, which is small enough, may undergo a quantum phase. 

Let us take a look at our Universe. For simplicity we will assume that the Universe 
has been filled with radiation from the big bang up to present. Since we are exploring 
the extension of the classical phase, we may use classical equation of motion. In 
addition, we set the intrinsic curvature to zero, as the current observations suggest it 
has not played a significant role in the evolution of the Universe so far. Then V9 = 

VoOo ■ — , where the subscript o refers to the present value. Let us fix Vq ~ 10 81 m 3 
to be the size of the observable Universe today and 6q ~ 10~ 11 c~ 1 s~ 1 . Let the unit 



of /3(v) be m and Ineq. (60) is saturated for 

(61) — « l(T 62 /3(z/) , 

ao 

hence the size of the Universe on the brink of the quantum era cannot be smaller 
than V ~ 10~ 186 /3 3 (^) m 3 « /3 3 (V)Zp, while the expansion rate is at least 9 w 
10 113 c _1 s _1 /3 _2 (^) ~ 10 58 tp 1 c _1 /3 _2 (^). Both these values may in principle be deep, 
deep into the Planck era. In particular we estimate that the energy density of matter 
was at that time: 

(62) p « 10 113 p„/r 4 (z/) • 

Now, we may expect that the quantum phase does not start prior to the Planck era 
and so /3(v) < 10 33 m 2 . Well, this is not strictly a constraint, because /3(v) can be 
by any order of magnitude larger provided that the Universe is respectively bigger 
than its observable part. Obviously, the above consideration relies on the simplest 
possible model of the big bang singularity. Nevertheless, it indicates that the Planck 
scale argument should be applied more carefully. As shown it is not applicable to 
the size or the energy density of the universe alone. One may however wonder if 
in another clock variable frame we would still obtain the same estimation for the 



domain of validity of classical dynamics. In other words, if inequality (60) would 



remain the same upon a different choice of clock, let us say - the size of the universe. 



At the moment we only notice that the equality (60) is only weakly sensitive to the 
value of w within its range. 



6.3. Modified Friedman equation. Putting the lower symbol of the Hamiltonian 
(54) back to the constraint equation (10) we obtain the effective Hamiltonian con- 
straint. In order to get the modified or semiclassical Friedman equation, we need to 
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divide the constraint by lapse iV = q 2w K l w ) and by volume V = q 1 - w . We getj^] 

3(1 -wfa 2 P K x {v) 2 + K (AK /a 

a w p g i-« + — — — j + 6A; g 

64RT {v)K 2 (u)q^ K G {v)K x {v) 



3(l-u>) 



-2 



1 + w 



p T q i-w . 

Now, we combine the definitions of (q,p) and (a,p) to arrive at: 

a\ 2 22/-, n2^(^) t-./ \kc 2 8nG 
a) +«l-f^ + iiW ? = ^, 

where 

1 ; 32^o(^)^2(^) 1 ^ K^K^v) 
As the result of quantization we obtain two corrections to the Friedman equation. 
First, the repulsive potential, which depends on the volume. We notice that as the 
singularity is approached a — > 0, the repulsive potential grows faster (~ a -6 ) than 
the density of fluid (~ a~ 3 ^ 1+w ^) and therefore at some point the contraction must 
come to a halt. Second, the curvature becomes dressed by factor B(v). This effect 
could in principle be observed far away from the quantum phase. However, we do not 
observe the intrinsic curvature neither in the geometry nor in the dynamics of space. 
Nevertheless, we note that B(v) ~ 1 for large enough v. The function A{y) 
for large v. Assuming v is large enough we obtain: 

d\ 2 22/ \2 v 1 87rG 



— 128 



(63) |-] +^(1-^ + = p, 



The form of the repulsive potential featuring in Eq. (63) does not depend on the 
state of fluid, which fills the universe. Therefore, we may conclude that the origin 
of singularity avoidance is quantum geometrical. Although the potential's coefficient 
weakly depends on the matter content through (1 — w) 2 , the dependence on w can 
be absorbed in the definition of v. The potential provides a kind of hardness to the 
collapsing space that resists its contraction. As the space contracts the gravitational 
interaction grows and increases the contraction rate even more. Then, at some point, 
the potential turns on and makes the contraction slow down until the space comes to 
a complete halt and rebounds. After the rebound the expansion initially accelerates 
and the potential turns off shortly after. Afterwards the dynamics becomes classical. 



5 The lapse and volume are functions of physical degrees of freedom and as such are considered 
as quantum operators. Thus, one might want to use lower symbols of N and V instead of classical 



counterparts to perform the division. However, we have already obtained semiclassical limit in (54 1 
and we should work with it by treating all the quantities, including N and V, classically. 
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6.4. Constraint from cosmography. Szydlowski et al. [27] have studied a cosmo- 
graphical constraint for the additional source term of the form fl mod (z + l) 6 to the 
Friedman equation of ACDM model. It corresponds to the geometrical correction 
in our modified Friedman equation. Interestingly, this kind of modification appears 
outside our construction in the framework of loop quantum cosmology, brane-world 
scenarios and others. The authors derived a constraint for the value of fi mo d- The 
constraint comes from a varied set of data: SNe la and radio galaxies sample, baryon 
oscillation peak (from Sloan Digital Sky Survey) and the so-called CMB shift pa- 
rameter. The result can be translated into our model's free parameter as: 



This constraint is, as expected, milder than the one we got from assuming that the 
correction comes into a play only after the Universe has entered the Planck era. The 
authors of j27J also estimate the constraint from the big bang nucleosynthesis to be: 



which is still a much weaker constraint then the first one that we have arrived at. 

6.5. Conclusion. The core result of this work is the occurrence of repulsive poten- 
tial providing a mechanism for the singularity resolution. The potential is quantum 
geometrical in nature and prevents the space from reaching the singular state. Even- 
tually, the space-time rebounds. The potential is a generic and unexpected feature 
of our quantization scheme. The coupling constant of the potential depends on the 
choice of fiducial vector, which is used to generate coherent states. We fully param- 
etrize this freedom with a positive and otherwise arbitrary parameter v. 

The free parameter is, in our mind, a desired result. First of all, v may be used to 
fit our model into observational constraints. We may also see the free parameter(s) 
as "modeling" the more fundamental structure of Quantum Gravity. On one hand, 
our theory is a merger of the principles of quantum mechanics and the dynamics 
of gravitational fields. On the other hand, unlike Weyl quantization, it includes 
the many ways in which one may do the merger. It encompasses the absence of 
observational clues of phenomena involving quantized gravity as well as the lack 
of knowledge of the fundamental principles on which the full Quantum Gravity is 
perhaps to be achieved one day. Our line of research may contribute to discovery of 
this structure. 

Our approach is bottom-up: we begun with quantizing homogenous and isotropic 
models. Next we will move to more complex ones. The natural extension of this 
work will be to allow for anisotropic evolution of space, which admits more complex 




from which we get: 




wfv = 2/3» < 10 131 m 4 



4(1 - w) 2 z/ = 2/3» < 10 120 m 4 
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singularity. We also plan to include different choices of clock variable and complement 
reduced phase space quantization with the Dirac method. 

Appendix A. Definition and essential properties of the affine group 

In this appendix, we leave aside all physical dimensions. Consider the affine trans- 
formation on the real line: 

R 3 t ->■ (a,b) -t := at + b 

where the pair dilation-translation parameters (a, b) G IR 4 ^ x KL 
The transformations (a, b) form a group with the composition rule: 

(a, b) (a , b') = (aa! , ab' + b) . 

The neutral element is (1,0) and the inverse of (a, b) is (a, b)' 1 = (1/a, —b/a). This 
group is called the affine group of the real line and denoted by the symbol Aff+(R). 
It has a left-action and a right-action on itself. The left-invariant measure is 

dadb 



whereas the right-invariant one is 

dadb 
a 

In the main text, we use a slightly different realization of the group, namely we 
perform a transformation of coordinates (a, b) — > (q, p) such that the measure is of 
the Lebesgue form, dqdp. One can easily check that this coordinate transformation 
is 

q = l/a, p = b 

for the left-action, and 

q = ln(a) , p = b 

for the right-action. This coordinate transformation affects the composition rule as 
well: 

(q, p) (q', p') = (qq 1 , ^ + p^j , (q, p) (q', p') = (q + q' , e q p' + p) , 
for the left-action and right-action respectively. 

In the main text we consider only the left-action. Also, the choice of coordinates q and 
p is such that they can be thought of as the configuration variable and its conjugated 
momentum parametrizing the phase space of FLRW cosmology, which is the half- 
plane II + = {(q,p) £ M 2 : q > 0}. In particular, since (q,p) G II + are canonical 
variables, the symplectic structure wrt these variables is of the diagonal form, and 
thus it follows that the measure is Lebesgue. That is why we did the transformation 
from the mathematically more natural (a, b) to the physically meaningful (q,p). 
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We already saw how Aff+(R) acts on K: in term of the new variables, 

s ->• (q,p) ■ s = - + p. 

q 

Let us make explicit Hilbert spaces, on which the (left-)action of Aff + (R) is unitarily 
represented. They are parametrized by a e Mfl 

H a := L 2 (R* + ,dx/x a+1 ) . 

The fact that these Hilbert spaces are based on the positive half-line, R+, is sug- 
gested by the usual situation in QM: the wavefunctions are maps from the classical 
configuration (in our case, the space where q takes values) to C. The requirement 
that such functions are square- integrable is also a standard QM requirement, while 
the weight in the Lebesgue measure is introduced for more generality. In the concrete 
study of physical operators, we conveniently choose a = —1, so that the measure 
becomes the usual one. 

The action of the operators U a (q,p) on ip 6 % a is defined as 

(U a (q,p)^)(x) = q a / 2 e i ^(x/q). 

It is easy to check that U a (q,p) is unitary. Its irreducibility has been shown in [H] 
(see also [T5]). 

In the main text, we use the action of U a (q ,p ) on a chosen ip to produce the so- 
called coherent state, \qo,Po), peaked on a classical phase space point (qo,Po)- As a 
function of x it is defined as 

(x\q , Po ) = (U a (q ,p )i/;)(x) = q a J 2 e^{x/q Q ) . 

Due to the affine group composition rule, we see that acting with U a (q,p) on \qo,po) 
produces 

(x\U a (q,p)\q ,p ) = (x\U a (q,p)U a (q ,p )\ip} = (x\U a (qq ,p /q+p)\ip) = (x\qq ,p /q+p) , 
i.e. 

u a (q,p)\q ,Po) = \qqo,po/q + p) =■ \q',p')- 

This shows that the coherent state \qo,Po) transforms covariantly under the action 

of unitary operator U a (q,p). 

At this point, Schur's lemma applies. 



It has to be said that the quantum theories arising from two Hilbert spaces of this family are 
unitarily equivalent, so one can effectively choose a at will, as it cannot play any role on the physics. 
In particular, in the main text the choice is a = —1, as this corresponds to the usual Schrodinger 
representation. 
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Schur's lemma. Let G be a group with U its UIR on a vector space V. If M is an 
operator on V such that U(g)MU(gY = M for all g G G, then M is a multiple of 
the identity 1 on V : M = c ■ 1. 

In our case, it is easy to check that the operator f dqdp\q,p){q,p\ satisfies the hy- 
potheses of the lemma. Therefore, it follows that 

J &qdp\q,p){q,p\ = c- 1 

Combining this operator identity with the projector \ip) (ip\ on a suitably selected unit 
vector and taking the trace allow to compute the constant. It is finite because 
the UIR is shown to be square integrable. This is nothing but the statement of the 
main text that the CS family {\q,p)} resolves the identity w.r.t. the measure dqdp/c, 
and is the starting point for coherent state quantization. 
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